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ON ; Abstract 



^ ■ A formalism is presented whicli allows covariant three-dimensional bound- 

state equations to be derived systematically from four-dimensional ones with- 
•/^ I out the use of delta- functions. The amplitude for the interaction of a bound 

state described by these equations with an electromagnetic probe is con- 
structed. This amplitude is shown to be gauge invariant if the formalism is 
truncated at the same coupling-constant order in both the interaction kernel 
CNJ , of the integral equation and the electromagnetic current operator. 
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I. INTRODUCTION 



The work of Fleischer and Tjon [|T]-^ over twenty years ago provided a foundation for 
use of the Bethe-Salpeter (BS) formalism for nucleon-nucleon (NN) scattering. Calcula- 
O I tions of electron-deuteron scattering using these BS wave functions were then performed by 
Zuilhof and Tjon Later work of van Faassen and Tjon |^,|^ developed a one-boson- 

exchange model for the coupled A^A^ — A^A system based on the BS formalism that provides 
a reasonable description of the experimental A^A^ phase shift data up to about 800 MeV. 
Given this progress in four-dimensional calculations, one might ask whether three- 
I dimensional calculations of the two-nucleon system are still useful. Reductions of the BS 
equation to a three-dimensional equation are motivated by several considerations. Firstly, 
three-dimensional calculations are simpler and hence allow the application of models for 
nucleon-nucleon interactions in three-body problems, e.g. the three-nucleon system and 
pion-deuteron scattering. In contrast, the only work using four-dimensional equations in 
three-hadron systems is an analysis for the three-nucleon system, by Rupp and Tjon 
which used separable interactions and a nonrelativistic treatment of spin degrees of freedom. 
An analysis of the three-nucleon system with the same ingredients as the BS analyses of the 
two-nucleon problem has not been realized due to the technical complexities involved. Sec- 
ondly, the meson-nucleon dynamics in A^A^ interaction models involves 'effective' degrees of 
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freedom that are designed for use in low-order truncations of the full BS kernel. If the degrees 
of freedom are truly "effective" then one cannot really regard the meson-nucleon coupling 
constants as fundamental (with the exception of the ttNN coupling constant). Given this 
point of view it is equally valid to construct an effective nucleon-nucleon potential in three 
dimensions, with the understanding that the coupling constants used will be different to 
those in the four-dimensional approach. Nevertheless, almost equivalent results for physical 
observables can be obtained in the three-dimensional and four-dimensional formalisms by 
making similar truncations of the scattering kernel and then adjusting the meson-nucleon 
coupling constants in both formalisms so that the empirical A^A^ phase shifts are well de- 
scribed. Thirdly, while it might appear that relativistic effects cannot be correctly treated 
in a three-dimensional approach, in fact all relativistic effects can, in principle, be studied 
in the three-dimensional framework. 

If relativistic effects are to be studied in a three-dimensional formalism some systematic 
reduction from the four-dimensional Bethe-Salpeter equation to a three-dimensional integral 
equation must be employed. Furthermore, if electromagnetic observables are of interest, this 
reduction should be implemented so that gauge invariance is realized in the three-dimensional 
approach. In this paper we present a three-dimensional formalism that is covariant and, after 
consistent truncation of the current operator and interaction kernel, is gauge invariant at 
any desired order. 

Consider the Bethe-Salpeter equation (BSE) for the four-dimensional, covariant, two-to- 
two amplitude T [|10|- |1^] , 



T = K + KGoT, (1.1) 

where Go is the free two-particle propagator, which in our convention is 

Goip[,P2■,Pl,P2)=^i2nfS^'\p[-Pl)6^'\p'2-P2)Goip■,P), (1.2) 

Go{p;P)=diipi)d2ip2), (1.3) 



with 



Pi - ei(Pi) p'> + ei{pi) - ir] 



where 



^iiPt) = ^Vi + mf, (1.5) 

for spin-zero particles. 



A^^(P^) \ ±;^fp\^y 7rP.+>n, ^ for spiu-half particlcs. ^^'^^ 
Here we have used the standard center-of-mass and relative four-momenta: 



Pi = PiP + p] P2 = T^2P-p, (1-7) 

with Ui + U2 = 1. Throughout this paper, when functions are written with P and p as their 
arguments they do not contain four-momentum conserving delta functions. In comparison, 
when written with pi and p2 as their arguments they do contain such delta functions. In 
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Eq. ( |1.4| ) ?7 is a positive infinitesimal. In Eq. ( |1 . 1|) , K is the Bethe-Salpeter kernel. In 
principle, K should include all two-particle irreducible two-to-two Feynman graphs. Solution 
of ( |1.1| ) with the full two-particle irreducible kernel is impractical and usually the kernel is is 
truncated to lowest order in the coupling constant, using the ladder approximation, K = V . 

A simple way to obtain an approximate three-dimensional equation from Eq. ( |1 . 1| ) is 
to ignore the relative-energy dependence of the interaction V in the loop integral of the 
ladder BSE. This leads directly to the Salpeter equation [|r^ , which, if only positive-energy 



states are considered, is equivalent to the equation of Blankenbecler-Sugar [|I6[ and Logunov- 
Tavkhelidze [0] 

T = Vinst + Mn.t(Go)r, (1.8) 

with the three-dimensional propagator 

(^0)= J^Go{p;P). (1.9) 

Tjon and collaborators have analyzed various hadronic systems, including elastic and in- 
elastic electron-deuteron scattering p8|-p0|, proton-proton bremsstrahlung ||2T[, and the 



quark structure of hadrons ||24-|24|, using a version of this equation with a modified three- 
dimensional propagator. These calculations have provided much insight into the role of 
relativity in few-body systems. 

In order to avoid the instant approximation in the three-dimensional formalism, it is 
necessary to correct for the difference between the full four-dimensional interaction V (or 
K) and the instant version thereof. Moreover, Eq. (|1.8| ) is usually obtained in the center-of- 
mass frame of the two-body system and the manner in which the amplitude is to be boosted 
to a frame other than the cm. frame must be defined. Much work to address these issues 
has been based upon the quasi-potential formalism in which a delta function constrains the 
relative energy of the interacting particles. Here, we discuss a different reduction to three 
dimensions that does not possess the unphysical singularities which arise in quasipotential 
reductions. 

A systematic formalism is found by splitting K into two pieces, one of which, Ki, does 
not depend on the zeroth component of relative four-momentum, and the other of which, 
K — Ki, does. Straightforward manipulations then lead to the following set of equations 
(see for details): 

r=[l-{K-K^)Gor'T,; (1.10) 
T, = K^{g)T^■, (1.11) 
g = Go + GoiK-K^)g; (1.12) 

{g)^J^-^gip'„p„p). (1.13) 

There is now a general way to account for all the relative-energy integrations that reside in 
{Q). Demanding that Ki be chosen such that Ki{Q) = Ki{Gq) produces 

ri = iri(Go)ri; (1.14) 

K, = {Goy\G^Kg){Go)-\ (1.15) 



3 



This choice of Ki is closely related to that used in the work of Klein [p^ -p9[] on three- 
dimensional reductions of four- dimensional equations, and to standard time-ordered pertur- 
bation theory. 

This formalism was developed in a non-covariant way in Ref. 



25| . In this paper we 



first show, in Section 0, how a simple modification of this argument allows one to obtain a 



covariant formalism. (Note that in contrast to Ref. ||2^ we do not discuss the approximate 
inclusion of crossed-ladder terms in the kernel of the integral equation here, although it 
is straightforward to apply the discussion of Section |I| to that aspect of Ref. also.) 
In Section |T| we take up the main issue of the paper: construction of a gauge-invariant 
amplitude in this formalism. We review the analysis of gauge invariance for the Green's 
function of the Bethe-Salpeter equation: 



G = Go + GqKG. 



1.16) 



Gauging such four- dimensional Green's functions is relatively straightforward. The critical 
issue is whether gauge invariance can be maintained within a consistent reduction scheme 
to three dimensions. We show how this is done in our approach to three-dimensional re- 
duction, thereby demonstrating that gauge invariance can be maintained order-by-order by 
consistently truncating the kernel of the bound-state equation and the bound-state current 
matrix element. 



II. COVARIANT BOUND-STATE EQUATIONS IN THREE DIMENSIONS 



In this section we generalize the procedure presented in Ref. ||2^ to yield a covariant 



procedure for obtaining three-dimensional bound-state equations. While the technique is 



essentially that discussed in Ref. we display it here for completeness and clarity. 

Consider the Bethe-Salpeter equation for the bound-state vertex function, which we write 

as: 

r{p;P) = 1 0^^Kip,p';P)Goip';P)Tip';P). (2.1) 

As discussed in the Introduction we now split the kernel K via 

Kip,p';P) = K,ip,p';P) + K2{p,p';P). (2.2) 

However, in contrast to the case discussed in the Introduction, here we make this decom- 
position in a covariant way. To do this we must write the relative four-momenta p and p' 

as 

p = P\\dD + P±d; p' = P\\dP> + p'±d- (2.3) 

Here, for any four-vector k, k\\ = k - D, with D a unit four-vector in the direction of D. The 
piece of the kernel Ki is then chosen so that it does not depend on p\\£) or p'l^. 

The choice of vector D is completely arbitrary and nothing depends upon it if the kernel 
used in the three-dimensional equation is treated exactly. Ultimately, of course, the kernel 
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must be truncated, and this introduces some dependence on the vector D into the three- 
dimensional formahsm. This will be discussed further at the end of this section. 
If the non-covariant vector 

D = (1,0,0,0), (2.4) 

is chosen, then the reduction to three dimensions is exactly that developed in Ref. [^. In 
order to obtain a covariant formalism for the bound-state problem, D = P must be chosen, 
because P is the only available four-vector of the problem that is conserved. With this 
choice, Eq. (|2.2|) becomes a manifestly covariant decomposition of the kernel K, since then 
it is expressed solely in terms of Lorentz-covariant objects. Obviously, this choice coincides 
with Eq. (p.4|) in the center-of-mass frame of the two-body system. But, in all other frames. 



this decomposition differs from that of Ref. ||25|, since it boosts from the cm. frame in a 
covariant way. We keep D distinct from P in what follows in order that other choices of D 
may be considered for current matrix elements. 
If Fi is now defined by 

Up;P) = I -0^,K,ip,p';P)Goip';P)T{p';P), (2.5) 

it follows that, since Ki does not depend on py^), neither does Fi. It is then a simple matter 
to show that Fi obeys the equation: 

ri(pxD; P)=l '^-0^K,{p^D,p%- P)g{p',p"; P)Fi(p'l^; P), (2.6) 

where the four-dimensional Green's function Q is defined by the integral equation 

g{p',p"; P) = Goip'; P)5^'\p' - p") + j J^Mp''^ P)K2{p' ,p"'; P)g{p"\p"; P). (2.7) 

Since Q is the only four-dimensional object appearing on the right-hand side of Eq. ( ^.6p , 
the equation reduces to one involving only three-dimensional integrals: 

ri(pXD;P) = / '^^^j0^K,{p^D,p'^^-,P)mP'^Dy^D^^ (2.8) 

Here, and in the rest of the paper, the notation {Q) denotes integration of the Green's func- 
tion Q over four-momentum components parallel to D, in order to yield a three-dimensional 
quantity: 

mp^D.p%; P)^j ^^j0^Gip',p"; P)- (2.9) 



In order to obtain F from Fi the equation 

F(p;P) = T,{p^n;P) + J ^^K2{p,p'-,P)G,{p'-P)T{p'-P). (2.10) 
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must be used. The T obtained from solving Eq. ( p.8| ) and using Eq. ( p.lO| ) is identically 
equal to that found by solving the original Bethe-Salpeter equation ( p.l|) . The difficulty 
is that the integral equations ( p.7| ) and ( p.lO| ) are themselves four-dimensional. So, at this 



stage we have not simplified the problem at all. 

We now consider whether there is a choice for Ki which simplifies the form of Q. In 
particular, we seek a Ki such that the propagator (Q) in Eq. (^78|) may be replaced by the 
free propagator (Go). In other words, the contribution of the last term of Eq. ( ^.7| ) must 
vanish once the integrals over pj^ and are performed. Writing K2 = K — Ki then leads 
to the defining condition: 

{Go){j>^d;P)K,{p^d,vW,P){Go){vW.P) = { J ^,Go{p;P)K{p,p";P)g{p",p';P)) 

(2.11) 



with Q still defined in terms of Ki by Eq. ( p.7| ). With this choice of Ki Eq. ( p.8| ) involves a 
simple propagator: 

Up±n;P) = I '^^K,{p^o,p'^D\P){G,){p%-.P)V,{p%-P), (2.12) 

and all the complexity is transferred to the interaction kernel Ki. We write the equation 
( p.l2| ) formally as: 

ri = iri(Go)ri. (2.13) 

In order to calculate the kernel Ki defined by Eq. (|2.11| ), a perturbative expansion of 
the four-dimensional kernel, must be made: 

00 

K = Y,K^^'\ (2.14) 

2 = 1 

Then, to second order in the coupling constant we see that Ki is given by 

irf) = {G,)-\G,K^'^G,){G,)-\ (2.15) 
while at fourth order in the coupling constant we have 

= (Go)-^((Goir(')Go) + {GoK^''>GoK^'^Go)yGo)-'-Ki'\Go)K['\ (2.16) 

Here the propagator (Gq) is 

{GoKp^n; P') = J^-^ p ^ ' (2.17) 

P\\D + 11]- ei- 62 P\\D -ir] + ei + 62 

where Af = Af (pj_|_^). These operators are defined as in Eq. (|1.6|), but with ej now given by 



^iiPiLo) = v"^? ~ Pi±D- other words, e's and A's are now defined as functions of 

not pj. These definitions then only coincide with the usual definitions of these quantities in 

the two-body center of mass frame. 
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The form ( |2.17| ) reveals an inconsistency in Eqs. ( p.l5| ) and ( p.l6|) . In a spin-half theory 
the inverse of (Gq) does not exist, since the propagator has no components in sectors where 
one particle is in a positive-energy state and the other is in a negative-energy state. However, 
if we only take and use matrix elements of equations such as ( p.l5|) and (|2.16|) in ++ and 

states the results given above are entirely correct. This caveat must be borne in mind 

when examining the equations presented below. The difficulty springs from demanding 
that {Q) = {Go) when the entity on the right-hand side of this equation does not allow 
propagation in H — or — h states. This problem can be completely avoided by use of the 
modified free propagator of Refs. ]TB|- pU| , p2| , p3| , p5| , pU| -|5^ , for which {Go)~^ is uniquely defined 
in all p-spin sectors. The results of this paper may be extended to that case. However, in 
order to focus our discussion on covariance and gauge invariance, we have chosen not to 
discuss this extension of the formalism here. 

If we consider a spin-zero or spin-half field theory with a Yukawa interaction, then the 
interaction k[^^ takes the form displayed in Ref. up to trivial kinematic replacements: 

{Go)ip'^D; p)k?\p'^^,p^d; P){Go)ip±D; P) = 




where e' 



ml 



Pi±D' 



A. 



The interaction Ki can be systematically improved by including the next terms in the 
coupling constant expansion of Eq. ( |2.11 ). In general, the Ki defined by Eq. ( 2.11 ) is the 
two-particle irreducible interaction, with two-particle irreducibility defined with respect to 
the propagator (Go)- 

As in the Klein approach (see Refs. p6|- p9|] ) the three-dimensional Green's function 
(Go + GqTGq), where T is the Bethe-Salpeter amplitude, may be expanded as a series of 
diagrams. From the arguments given so far in this section, it is now clear that the Klein rules 
discussed in Refs. P5|-p7[], will only be modified by the following kinematical replacements: 



E 
Pi 

uj{p - p ) 



PiLD'i 

^iiPi±D) = \]mi-Pi%-. 

Uj{p±D - p'±d) = \//^^ - iP±D-p'±Dy 



(2.19) 

(2.20) 
(2.21) 
(2.22) 



7 



Any choice for Ki may therefore be expressed as a set of diagrams and calculated in a 
straightforward way. Once again, we observe that if the non-covariant choice (|2.4| ) is made 
for then these replacements have no effect. Since the Klein rules coincide with those of 
time-ordered perturbation theory in ++ states, an interaction generated using a particular 
set of diagrams in the formalism discussed here with this choice for D will agree with the 
time-ordered perturbation theory result for the same diagrams, provided only positive-energy 
states are considered. 

However, unlike either the time-ordered perturbation theory or Klein interactions, Ki 
may be defined so that it boosts in a covariant way. Consider, for instance, the choice 
D = P. The result of the above procedure is then a covariant three-dimensional interaction 
for use in Eq. ( p.l2|) . In the cm. frame, where P = (1,0,0,0) this interaction agrees with 



the Klein result. After solving for Fi in the cm. frame we can simply make the replacement 
p p±p and thereby obtain a covariant vertex function in any frame. Once this equation 
is defined in one frame a kinematic boost defines it in all other frames. 

If D is not chosen to be parallel to the total four-momentum of the two-body system 
then vertex functions in an arbitrary frame cannot be obtained from a kinematical boost of 
the rest-frame vertex function with the same D. To see this, consider some arbitrary choice 
for the four- vector D. Suppose that P is chosen to be a four- vector describing the deuteron 
at rest and that we wish to calculate Ti{k±D', P')- Under a Lorentz transformation C 

T,{[Ck]Dr, CP) = S{C)r^{k^n; P), (2.23) 

where D' = CD, S{C) = Si{C)S2{C), and Si is the unitary representation of the boost C 
acting on the spin-half particle i. Consequently, Ti{kj_£); P') can be related by a Lorentz 
boost to Ti{kj_Do; P), where C is chosen such that P' = CP, and Dq = C^^D. Therefore 
r can always be related to a rest-frame vertex function, but, unless D is chosen to be the 
total four-momentum of the two-body system, this rest-frame vertex function corresponds 
to a reduction of the four-dimensional Bethe-Salpeter vertex different from that used in 
calculating Tiik^j:,; P). 

When a truncated kernel Ki is used this difference matters, since then for each choice of 
D the kernel of the bound-state equation (p. 12 ) is different. Consequently, the bound-state 



equation must be re-solved at each value of P' in order to determine the vertex function 
r(^±D; P')- If a truncation of the kernel Ki is employed, the solutions for values of P' ^ P 
will not correspond to the same bound-state mass as the vertex function Ti{k±D; P). It is 
possible to minimize this effect by using a sufficiently high-order truncation of the kernel. 
A simpler expedient is to use a low-order truncation and correct this by rescaling the kernel 
as a function of D so as to maintain the same invariant mass for the initial and final states 
of Eq. ( p. 141 ) . Because of Eq. ( |2.2^ ) this rescaling may equivalently be viewed as a function 



of the total four-momentum P'. 



III. CONSTRUCTION OF A THREE-DIMENSIONAL GAUGE INVARIANT 

AMPLITUDE 

In Section || we showed how to derive a three-dimensional two-body bound-state equation 
from the corresponding four-dimensional Bethe-Salpeter equation. With the choice D = P 
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this method maintains the Lorentz invariance of the original theory. The calculations of 



Refs. 1 25. 33] indicate that this formalism holds promise for accurately approximating the 
results of Bethe-Salpeter equations. However, in any reduction of a gauge- invariant, four- 
dimensional theory to three dimensions, a critical test of the reduction procedure is whether 
it can consistently maintain the gauge invariance of the full four-dimensional theory. In 
this section we first consider the Green's function for the interaction of two particles with 
the photon constructed from the Bethe-Salpeter Green's function. The general construction 
of currents guarantees that this four-dimensional Green's function obeys a Ward-Takahashi 
identity (WTI), as required for gauge invariance. (Obtaining such a WTI in the four- 
dimensional theory in general requires the inclusion of meson-exchange currents.) We then 
turn our attention to the corresponding three-dimensional formalism, and show that truncat- 
ing the two-body current and the two-body interaction Ki at the same order in the coupling 
constant leads to a formalism for the calculation of electromagnetic processes which is both 
gauge invariant and Lorentz invariant. 

We begin our discussion by reviewing the WTI of the Bethe-Salpeter Green's function. 
Consider the two-body Green's function G that is the solution of the Bethe-Salpeter equation 

G = Gq + GoKG. (3.1) 

Define the Green's function for the interaction of a free single particle i with a photon of 
fixed momentum Q via: 

^Sfe, Q) = dM + Q)3fmdi{Vi)- (3.2) 

Here the operator j^*^ is assumed to be such that this one-body photon Green's function 
obeys a WTI |3|,|§ 



Q) = e^{dM - d,{pi + Q)). (3.3) 
Defining the two-body analog of Eq. 

GoiPi + Q,P2)jl!\Q')d^\p2)Goipi,P2) + Go{pi,P2 + Q)d^\pi)j;-\Q')Goipi,P2) 

= GoiPi + Q:P2)jl!\Q')d^\p2)GoiPuP2) + (1 - 2) = G^ ^{pi , P2 , Q) ■ (3.4) 

Here, and throughout the rest of the paper, the notation (1^2) indicates that the momenta 
of the two particles must be swapped, and the labels exchanged. Therefore, the (1 2) 
pieces of any expression represent the photon coupling to whichever particle it did not couple 
to in the first part of the expression. We observe that the Green's function GqJ^ also obeys 
a WTI 

Q>'G^o^{p^,P2,Q) = ei(Go(pi,P2) - Goipi + Q,P2)) + {1 ^ 2). (3.5) 

Now let be the Green's function for the interaction of one photon with the interacting 
two-particle system. Note that in in the two-nucleon system the quantities Cj will include 
isospin operators, so care must be exercised in ordering charges and any interactions which 
also involve isospin operators. We may write the following equation for a gauge invariant G^, 
by simply allowing the photon to be inserted anywhere on the right-hand side of Eq. ( p.l|) . 
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G 



, Gl^ + Gl^KG + G.KGl + G.KIG. 



(3.6) 



Here K'J^ is found by coupling the photon to every internal charged line in the kernel K. 
Using Eq. ( |3.1| ) and the definition ( |3.4|) , Eq. ( p.6| ) may easily be solved for 

Gl{k[,k',-Kh-Q) = J ^-^^G{k[,k',;p, + Q,p2)/^HQ')d^\p2)Gip^,p2-,h,h) 



+(1 ^ 



2) + / 



(2>r) 

d'^p'i d^p'2 d^pi d'^p2 



{2n) 



16 



G{k[,k'2;p[,p'2)KVp[,P2;Pi,P2;Q)G{pi,p2;ki,k2). (3.7) 



(Note that there is an overall delta function 6^^\k[ + k'2 — ki — k2 — Q) on both sides of this 
equation.) Using the identity ( |3.5|) in Eq. ( p.7|) and the explicit form of Eq. (|0| ) we find 



Q^Gl{k[, k'2; k,, k2; Q) = e^G{k[ - Q, k'^; k,, k2) - G{k[, k'^; k, + Q, k2)e, 
d'^p[d'^p'2d'^Pid'^P2 



+ 



+ 



(27r)i6 
d'^p\d'^p'2d'^Pid^P2 



{2n) 



16 



G{k[,k2;p[,P2) [K{p[,p2;pi + Q,P2)ei 

-eiK{p[ - Q,p2;pi,p2)]G{pi,p2;ki,k2) + (1 ^ 
G{k[,k'2;p[,P2)Q^K^{p'i,P2]Pi,P2]Q)G{pi,p2;kuk2). 



(3.8) 



This reduces to a WTI for G^, 



Q''Gl{k[, k'2; ki, k2] Q) = e^G{k[ - g, k'2] ki, k2) - G{k[, k'2; ki + Q, k2)ei + (1 ^ 2), (3.9) 
provided that 

eiK{p[-Q,p2;pi,p2) - K{p[,p2;pi + Q,p2)ei + {1 ^2) = Q^K'l{p[,p2;puP2;Q), (3.10) 

which is the WTI for the interaction current. (Similar identities are used in the construction 
of a gauge invariant electromagnetic matrix element for the Gross — or spectator — formalism 
in Refs. p6|j37[1 .) The result ( |3.10| ) is completely general, and will always hold if the two- 
body current K'^ is constructed in a gauge invariant way. In Appendix we show how 
Eq. (p.lO|) is achieved in the special case where the photon couples only to particle one, 
i.e. the exchanged particles and particle two carry no charge. By contrast, if the particles 
are, for instance, nucleons carrying isospin which exchange isovector mesons, the analysis is 
more complicated, partly because the isospin operators in K and G do not commute with 
those in the charges ei and 62- However, as alluded to above, the additional terms generated 
because of this non-commutativity are canceled by terms arising from the coupling of the 
photon to isovector mesons. (The interested reader may consult Ref. for details on this 
problem.) 

Now, in the case of a bound state we have the following decomposition of the Green's 
function 



G{k[, k'2, ki, ^2) = 

5^'^{k[ + k'2-k^-k2) 



GoiK, k'2)^^^^^0^^^^Go{k„ k2) + R{k[, k'2, fci, k2) 



(3.11) 
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where P = ki + k2 = k[ + k!^ is the conserved total four-momentum of the two-body system, 
M is the two-body bound-state mass, and the piece R is regular at = M^. We then insert 
Eq. ( |3.11| ) into ( p.7| ) and extract the double-pole part of the resulting expression, in order to 
obtain the amplitude for interaction of the bound-state with a photon of momentum Q |p8|| . 
Expressing the result in terms of total and relative four-momenta yields: 



A,{P,Q)= J -^^T{p + Q/2;P + Q)G2^{p;P,Q)T{p;P 



+ J ^-^r{p';P + Q)Go{p'-,P + Q)K;{p',P + Q;p,P-,Q)Go{p-,Pnp;^ (3.12) 

From Eq. (|3.5|), Eq. (|3.1CI| ), and the bound-state BSE, Eq. ( p.l|) , it is straightforward to show 
that 

Q^A,{P,Q) = 0, (3.13) 

as required by gauge invariance. 

In the context of the formalism of Section || the question is how to maintain this gauge 
invariance when the reduction to three dimensions is made. In particular, in Section ^ we 
showed how Fi, which obeys a three-dimensional integral equation, could be used to find F. 
Inserting the expression ( |2.1CI| ) into Eq. ( p.l2| ) and using the definition of Q we find 

d^k'd^p'd^pd^k- 



A,{P, Q)=J ^ly, Tiik'^D-, P + Q)<5{k\p'- P + Q) 



i2ny 

' j^;\Q')d',\v,P - p)6('\p' -p- Q/2) + j(^\Q')d^\u,P - p)S('\p' -p + Q/2) 
+K;{P + Q,p'- P,p- Q)] g{p, k- P)T,{k^n; P). (3.14) 

An important part of Eq. ( |3.14| ) is that we use the same D in the initial and final states. 
From the discussion of Section it might appear that the natural choices are D = P in 
the initial state, and D = P + Q in the final state. However, if these choices are made 
then the formula ( |3.14 ) is actually slightly misleading, because the two vertex functions so 



defined depend on pieces of the four-vectors k and k' which reside in different subspaces. 
Thus, for instance, the two ^'s appearing in Eq. ( p. 141) are actually different objects. This 
immediately leads to significant complications. In particular it makes it very difficult to 
define a three-dimensional gauge invariant approximation to the amplitude ( p.l2| ). 
In this current matrix element, the choice 

D = {P + P')/2, (3.15) 

may be made in both the initial and final state. With this choice the vertex functions 
^i{k±D,P + Q) and Ti{kj_u, P) cannot both be obtained from a kinematical boost of the 
rest frame wave functions with D = (1, 0, 0, 0). However, in spite of the resultant variation in 
the bound-state mass discussed in Section ||, the matrix element ( p.l4 ) is covariant provided 



that both the vectors P and P' are boosted using the same Lorentz transformation C This 
is easily shown by expanding Eq. ( |3.14|) in powers of K — Ki and then using Eq. ( p.23| ) , the 
definition of Ki, and the known covariance properties of K and Gq. 
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The choice ( p.l5| ) is a useful one because it allows the straightforward definition of a three- 
dimensional matrix element for the interaction of the bound state with an electromagnetic 
probe which is not only covariant, but is also gauge invariant. To show this, first note that 
in order to solve the integral equation (|2.8| ), we in fact choose Ki so as to impose 



{G) = (Go). (3.16) 

If this condition is imposed to all orders in the coupling constant then the result is an infinite 
series for Ki . If we impose the condition ( ^.16] ) order- by-order in the expansion in — Ki , 
the condition defines Ki to the same order. Truncation of the kernel is necessary for a 
practical analysis and it becomes crucial to inquire if a corresponding approximation for the 
matrix element ( |3.14| ) exists that maintains the gauge invariance of the theory. It turns out 
that the current matrix element ( p.l4|) is gauge invariant if Q{J^ + K'^)Q on the right-hand 



side of Eq. ( p.l4| ) is expanded to a given order in the coupling constant and the kernel Ki 
used to define Fi is obtained from Eq. ( p.l6| ) by truncation at the same order in the coupling 
constant. 

To see this we first abbreviate Eq. ( |3.14| ) as: 

A,{p, Q) = f i(p + Q)g{p + g)( + Ki)g{p)v^{p), (3.17) 

where all integrations over relative momenta are now suppressed. We split this into two 

pieces, one due to the one-body current J^, and one due to the two-body current KJ^. 

(2) 

Suppose that Ki has been truncated at lowest order, i.e. Ki = . Then, in the piece, 
we expand the ^s and retain terms up to first order in K — Ki. To this we add a piece from 
the two-body current, in which we stop this expansion of Q at zeroth order in — Ki, i.e., 
write Q = Gq. Thus, we define a first-order approximation to A^, A^^\ by 

= f i(P + Q){Gl^)T,{P) + f i(P + Q){Go{P + Q){K{P + Q)-K,{P + Q))Gi;)V,{P) 
+ f i(P + Q){Gl^{K{P) - K,{P))Go{P))UP) 

+ T,{P + Q)(Go(P + Q)K;Go{P))T,{P). (3.18) 

It follows from Eq. ( ^1231) that is a covariant matrix element. To show A^^^ is also gauge 
invariant we contract with the four-vector Q and use the identities p.5| ) and ( |3.1CI| ) This 
leads to 

QM« = ri(P + g)ei(Go(P))ri(P) - f i(P + g)(Go(P + Q))e^T^{P) 

+ f i(P + Q){GoiP + Q){K{P + Q)-K,{P + Q))e^GoiP))T^{P) 

- f i(P + Q){GoiP + Q){K{P + Q)-K,{P + Q))Go{P + Q))eiT,{P) 
+ ri(P + Q)ei(Go(P)(ir(P) - K,iP))Go{P))T,iP) 

- f i(P + g)(Go(P + Q)ei{KiP) - K^iP))GoiP))T^iP) 

+ f i(P + Q)(Go(P + Q){eiK{P) - K{P + Q)ei)Go{P))T,{P) + (1 ^ 2). (3.19) 

The bound-state equation (|2.12|) can now be used to show that the terms in the first two 
lines cancel the terms in the last two lines, and so: 

gM« = -f i(P + Q)(Go(P + Q){K{P + Q)-K,{P + Q))Go(P + Q))elT^{P) 

+ri(P + Q)e,{GoiP)iKiP) - K,iP))GoiP))T,iP) + (1 ^ 2). (3.20) 
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But, at second-order in the coupling constant, Eq. (|3.16| ) leads to Eq. ( p.l5|) for Ki. It 



immediately follows that if Ki has been defined in this way, the corresponding amplitude 
for electromagnetic interactions of the bound state, as defined by Eq. ( p.l8| ), obeys 

g^!,') = 0. (3.21) 



It is straightforward to check that the same result holds if Eq. ( 2.12|) is truncated at 



next order, while the one-body current pieces of Eq. ( p.l4|) are expanded to second order in 



K — Ki and the two-body current pieces are expanded to first order inK — Ki. That is, if the 
interaction k'i'' defined by ( |2.16| ) is used, then defining a vertex function A^^'^ by expanding 
Eq. (|3.17|) to fourth order in the coupling constant leads to Q'^A^^'> = also. Thus, we 
conclude that truncating the kernel Ki defined by Eq. (|3.16|) and the electromagnetic vertex 



defined by Eq. (|3.17|) at consistent order in the coupling constant yields a gauge invariant. 



covariant, electromagnetic matrix element. 

It should be pointed out that A^^^ includes contributions from diagrams where the photon 
couples to particles one and two while exchanged quanta are "in-flight". These contribu- 
tions are of two kinds. Firstly, if the four-dimensional kernel K is dependent on the total 
momentum then gauge invariance requires the presence of terms representing the coupling 
of the photon to internal lines in K. Secondly, even if the kernel does not depend on the 
total four-momentum, e.g. it is a OBE kernel, terms arise in the three-dimensional for- 
malism where the photon couples to particles one and two while an exchanged particle is 
"in-fiight". These must be included if our approach is to contain a WTI. (See Fig. |I| for a 
diagrammatic interpretation of one such term.) Many years ago, in the course of work on 
deuteron photo-disintegration, Pearlstein and Klein ||39| showed that contributions such as 
these had to be included if Siegert's theorem (which is a consequence of gauge invariance) 
was to hold in the Klein approach. Our systematic derivation shows the general form of 
the additional pieces of the current operator which are required if gauge invariance is to be 
maintained. 

Gross and Riska and Kvinikhidze and Blankleider have discussed the construc- 



tion of gauge invariant amplitudes for the interaction of a photon with the two-body system 



within the Gross (or spectator) quasipotential approach ||0,|T|. It would be interesting to 
compare the systematic approach to constructing electromagnetic matrix elements developed 
here with the discussions of Refs. p6|j37| . 



IV. CONCLUSION 

In this paper, a systematic formalism is presented for the reduction of the Bethe-Salpeter 
equation to a three-dimensional form in which components of the relative momentum parallel 
to a direction D are integrated out. By choosing D = P, where P is the total four-momentum 
of the two-body system, a covariant formalism is obtained for a bound state. 



This part of the development follows closely the analysis of Ref. p5|, where the corre- 
sponding formalism based on the choice D = (1, 0, 0, 0) was constructed. However, here 
we simplify the analysis by omitting the treatment of crossed graphs in order to focus on 
covariance and gauge invariance. 

Consideration of the matrix elements of a conserved current shows that the gauge in- 
variance of the underlying four- dimensional theory is preserved in the three-dimensional 
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reduction. Moreover, it is possible to truncate both the infinite series expansion for the 
interaction kernel used to define the bound states and the expansion of the currents at the 
same order in the coupling constant, and maintain this gauge invariance. 

The three-dimensional reduction for current matrix elements is based upon the use of 
the same direction D in initial and final states. Here, the choice D = |(P + P') is suggested, 
where P {P') is the initial (final) total four-momentum of the two-body system. Other 
possible choices, e.g. D = P or D = P', could be used just as well, and may be preferable 
in some circumstances. The current matrix elements are covariant for any choice of D that 
is a linear combination of the two independent four-vectors P and P'. 

When truncations of the interaction kernel are employed in order to implement a practical 
calculation the results depend on the choice of D. The dependence on D in a truncated 
kernel, Ki, can be reduced systematically by keeping higher-order terms. Alternatively, the 
invariance of the mass of initial and final bound states can be enforced by a rescaling of the 
coupling constants in Ki. 

There are two important advantages of the reduction formalism based on integrating out 
momentum components presented here. Firstly, it is systematic. Secondly, as discussed in 
Ref. p5[, it avoids unphysical singularities. Such singularities generally arise in systematic 
quasipotential reductions to three dimensions. Typically such reductions involve splitting 
the propagator into a part involving a delta function and a remainder. Corrections to the 
theory have unphysical singularities, such as particle production thresholds at energies below 
the physical values. Moreover, current matrix elements obtained using quasipotential reduc- 
tions generally involve similar unphysical singularities which must be either disregarded 
or suppressed, with uncertain consequences. Such difficulties have motivated the present 
approach. This systematic formalism in which gauge invariance can be maintained by con- 
sistent approximations in the current operator and interaction kernel should prove to be a 
useful tool in the calculation of relativistic bound state properties. 
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APPENDIX A: PROOF OF WARD-TAKAHASHI IDENTITY FOR KJi IN A 

SPECIAL CASE 

In what follows we assume that the photon couples only to particle one. If the interaction 
between particles one and two is mediated by quanta which carry charge, then the discussion 
given here must be modified to incorporate terms where the photons couple to the exchanged 
quanta. For instance, if particles one and two are nucleons carrying isospin which exchange 
isovector mesons, the analysis is more complicated, due to the non-commutativity of isospin 
operators in K and G with those in the charges ei and 62- However, the additional terms 
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generated via this complication are canceled by terms arising from the direct coupling of the 



photon to the isovector mesons |36 



Consider the nth order piece of the kernel K. If 62 = and the exchanged particles do 
not couple to the photon we need only concern ourselves with the particle one pieces of the 
three terms in Eq. ( p.lO[ ), since the particle two and meson propagators may be arranged 
so as to take exactly the same form in all three pieces. Suppose then that the momenta on 
internal lines of particle one are labeled p^i\pf'\ ■ 



Pi^ to the right of the photon insertion, 



and Pi^ + QiPi^^^ + Q, • • • -\- Q to the left of the photon insertion. Meanwhile, the 

vertices governing the interaction of the exchanged quanta with particle one are labeled 



(1) 



Vi'^\ With these conventions we may write the particle one pieces of the three 



terms in Eq. ( |3.1CI| ) as: 



n-l 

K{p',.p'2.vi + Q.P2) : vt^ n + Q)yi'^ 



K{p[ - Q,p',;Pi,P2) : V;("^ n di{p':^)Vl 



i=l 
n-l 



i=l 



n-l 



i^;(p'i,P2;Pi>P2;Q):EW 



in) 



n— 1 

n d,{pf^+Q)Vl 
i=j+l 



9}>'j{p['\QM 



I[d,{p?)V^ 



i=l 



(Al) 
(A2) 
(A3) 



where the product symbols are defined to give one if the upper index of the product lies 
below the lower index. Note that the changes in particle one's external momenta guarantee 
that the particle two and meson pieces of the diagram will be identical in all three cases. 

Now calculating Q'^K'^ using the identity ( p.3|) we obtain, for the particle one pieces of 
this object: 



n-l 



(n) 



n-l 

n d,{pf + Q)Vl 
i=j+l 



eMp?)vl'^ 



n-l 



(n) 



n-l 

n d,{p? + Q)Vl 
i=j+l 



Udi{p?)v}'^ 

i=l 

di{p? + Q)eiVl'^ 



[ 

i=l 



i[di{p?)vr^ 



(A4) 



This is a telescoping series, as the jth term from the first piece cancels the j + 1th term 
from the second piece. The only structures which survive this cancelation give precisely the 
particle one piece of 



eiK{p[ - Q,p'2;pi,P2) - K{p[,p'2;pi + Q,P2)ei. 



(A5) 



Note that the series can only be summed in this way if the vertices commute with the 
charge operator ci. However, if the particles are distinguishable and [V^/"''',ei] = then this 
argument suffices to prove Eq. ( p.lO|) . 
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FIGURES 



FIG. 1. One example of a graph which is required in our formahsm in order to maintain gauge 
invar iance. 
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